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Abstract 

A relativistic model for the quasifree photoproduction of ry meson 
from complex nuclei is developed. The interactions between fields are 
introduced through effective Lagrangians. Contributions from several 
nucleon resonances as well as nucleon Born terms and vector meson 
exchange diagrams are included. Nucleon and rj wave functions are so- 
lutions of Dirac and Klein-Gordon equations, respectively. Final state 
interactions of the outgoing particles are included via optical poten- 
tials. The effects of these interactions on the cross sections and photon 
asymmetries are studied and are found to be large. Calculations indi- 
cate that at energies near threshold the exclusive reaction takes place 
mainly through formation of the Sii( 1535 ) resonance. Comparisons 
with the non-relativistic calculations show differences between the two 
approaches both for the cross sections and photon asymmetries. We 
give some detailed predictions for the reaction observables for exclu- 
sive photoproduction on ^^C, ^^O and ^''Ca. 
PACS number(s) 25.20.Lj, 24.10.Jv, 24.70-hs, 13.60.-r, 13.60.Le 
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1 Introduction 



r] meson production reactions explore a rich domain of phenomena at the 
interface of nuclear and particle physics. These reactions have the potential 
for expanding our understanding of the formation of nucleon resonances and 
their propagation in the nuclear medium. Moreover, as has been shown re- 
cently, these reactions have allowed more precise determination of the mass 
of the T) meson and the study of its rare decays Q . The spin-isospin char- 
acteristics of the 7] lead to a selectivity of resonance channels that can be 
formed through its interactions with a nucleon. At energies near the thresh- 
old of its production, there is a significant preference for the excitation of 
the Sii(1535) resonance which is known to decay nearly 50% of the time to 
an 1] meson and a nucleon. By contrast, the Sii(1650), which has identi- 
cal characteristics, has only a decay rate below 2%. The understanding of 
these dissimilar branching ratios is still an open question and may be a re- 
flection of important and subtle differences in the detailed substructures of 
these resonances. The selectivity mentioned above makes the t] meson pro- 
duction reactions an important tool for studying the physics of the Sii(1535) 
resonance. 

Until recently the production of t] mesons using photon beams has been 
largely confined to proton targets. The advent of high-duty-cycle electron 
accelerators such as those at the Jefferson Lab, Mainz, Bates and other lab- 
oratories, has opened up novel possibilities of performing production experi- 
ments on complex nuclei f^. 

The elementary reaction 7 + p — *■ rj +p has been the subject of several 
studies. Benmerrouche et al. 0] developed an effective Lagrangian approach 
to study this reaction. Contributions from nucleon resonances, t-channel 
vector meson as well as the nucleon Born diagrams are included. Calcula- 
tions were performed at the tree level and included eight free parameters. 
The S'ii(1535) was found to be the dominant contributor to the reaction at 
energies close to threshold. In addition, from comparisons with the data, the 
authors concluded that there is no clear preference for using pseudoscalar 
( PS ) or pseudovector ( PV ) coupling forms for the interaction vertices of 
the 7] meson with the nucleon and nucleon resonances. Bennhold and Tanabe 

1^ developed a somewhat different approach to r] photoproduction on the 
nucleon. They used a coupled channel isobar model, in which the reaction 
is related to ( 7, vr ), ( vr, 77 ) and { f], r] ) reactions. They included con- 
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tributions from three resonances: Sii(1535), Pii(1440) and Di3(1520). The 
data on the elementary reaction appeared to favor the PS form of the ?7NN 
vertex. The authors employed the resulting elementary amplitude to study 
coherent and incoherent photoproduction of t] mesons on nuclei. This model 
for the elementary process was later extended by Tiator et al. to include 
Born as well as p and uj exchange diagrams through effective Lagrangians. 
The resulting elementary amplitudes were subsequently used to calculate the 
PWIA (Plane Wave Impulse Approximation) amplitude of the photoproduc- 
tion reaction on light nuclei. 

The elementary amplitude mentioned above was also used by Lee et al. 
[|] to calculate the amplitude for quasifree photoproduction of rj mesons 
on heavier nuclei. They performed exclusive as well as inclusive reaction 
calculations using a nonrelativistic approach. The initial bound nucleon is 
described by a harmonic oscillator wavefunction and final state interactions of 
the outgoing particles with the recoil nucleus are introduced through optical 
potentials. Two different optical potentials were used for the rj mesons. The 
authors found that the photon asymmetry is insensitive to the final state 
interactions of the outgoing particles as well as the size of the target nucleus. 

In this paper we develop a relativistic model for the (7, rjp) reaction on 
nuclei leading to specific final states of the residual nuclei ( exclusive reac- 
tions ). The study of the photoproduction reactions on nuclei compliments 
the ongoing studies of the production on free nucleons. In addition to further 
investigating the reaction mechanism, studies on nuclei can help limit any 
existing ambiguities in the elementary process. Moreover since the produc- 
tion process is mediated by certain nucleon resonances and vector mesons, 
these studies open the possibility of investigating medium modifications of 
the properties of these hadrons in nuclei. As mentioned earlier, the avail- 
ability of a new generation of accelerators capable of producing high quality 
beams of electrons and photons, has led to increased interest in eta photo- 
production physics. Many experiments that were not possible in the past are 
now within reach. 

The present work is devoted to quasifree photoproduction. This type of 
reaction has the advantage of allowing measurements at small momentum 
transfer to the residual nucleus. This, in addition to enhancing the possibil- 
ity of increased cross sections, compliments the information obtained from 
coherent and incoherent { 'y, ">] ) reactions which by necessity involve large 
momentum transfers. 
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In the present study we follow a relativistic approach. This is motivated 
by the impressive successes of the relativistic mean field approach to nu- 
clear dynamics. We have recently studied the difference between relativistic 
and non-relativistic treatments in photonuclear knockout and quasielastic 
scattering of electrons [|T0|, and have found that there are important differ- 
ences between the two approaches. In the course of the present study we 
shall compare our results to those of the work of Lee et al. [§| . 

The main ingredients of the present treatment are as follows. The effective 
Lagrangian of Benmerrouche et al. is used for the interactions between 
fields. Contributions from nucleon resonances, t-channel vector mesons as 
well as the nucleon Born diagrams are included. In addition to the Sii( 1535 
) resonance, we also include three spin-| resonances ( Pii(1440), Sii(1650), 
Pii(1670) ) and one spin-| resonance ( Di3(1520) ). The nuclear wavefunc- 
tions are solutions of the Dirac equation with strong scalar and vector po- 



tentials in the spirit of the relativistic mean field theory of Walecka []TT], [12 
Calculations are carried out in the PWIA limit and also in the Distorted 
Wave Approximation (DWA) in which the final state interactions are taken 
into account. Bennhold and Wright have developed a similar model for the 



photoproduction of kaons on nuclei |13]. A short note on the present work 
was reported earlier [|14|. 

In the following section we derive the amplitude of the A( 7, ?7p )A— 1 re- 
action and calculate the relevant observables. The results of our calculations 
are presented in section III and our conclusions are given in section IV. 



2 Formalism 

The model for the ( 7, ?7P ) reaction on nuclei has the incident photon in- 
teracting with a bound proton through some process denoted by the blob in 
Fig.l. A proton and an rj meson are produced. The diagrams contributing 
to the blob of Fig.l, at the tree level, are shown in Fig.2. The Born dia- 
grams ( s- and u-channels of the nucleon propagator ) are shown in Fig. (2a) 
and Fig. (2b), the t-channel vector meson diagram is shown in Fig. (2c) and 
Figs. (2d) and (2e) are the nucleon resonance poles. 

The starting point in the present approach is a relativistic interaction 
Lagrangian for a system of photons, nucleons and mesons from which one 
obtains the transition amplitude for the A( 7, ?7p )A— 1 reaction. The am- 
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plitude then is used to calculate the observables for the reaction. 



2.1 The Interaction Lagrangian 

In the photoproduction of 1] mesons from complex nuclei, the reaction is 
taking place within the nuclear medium. The dynamics of the nucleons 
within the nuclear matter are described through the relativistic mean field 



Lagrangian of Walecka [0, |T^. The rj meson is described by solutions to 
the Klein-Gordon equation. The interactions of the fields involved in the 
production reaction depicted in Fig. 2, are described through the following 
interaction Lagrangian 

In the present work we adopt the pseudoscalar form for the t] meson-nucleon 
states coupling. This choice is based on the work of Bennhold et al. ^ and 
Lee et al. [|| which showed preference for this form in the production on 
protons. The interaction Lagrangians which involve the nucleon and spin-i 
resonances are cast in the following forms [Q: 

ll^nNN = - igni^i^i^v (2) 
C,NN= - ei^j.A'^i; - ^i^a'-'^^F,,, (3) 

CvNN= - g.i'i.W^-^i^cj.^i^v^'', (4) 

Cv,, = P^e^^x^F'^-'V^'^T], (5) 
'CrjNR = - ig^NRipTRf] + h.c. (6) 
£^^= - -RT'"'^F.^ + h.c., (7) 
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where 

for Sii resonances F = 1 and T^'^ = 'y^cr^'^ 

for Pii resonances F = 75 and T'^'^ = a^^ (8) 

and are the anomalous magnetic moment of the proton and the nucleon 
( proton ) resonance, respectively and M/j is the resonance mass. The tensor 
V^y{x) is related to the vector meson field by 

V^,{x) = d^V,{x) - d,V^{x) (9) 

The interaction of photons and 77 mesons with the spin-| ( Di3(1520) ) reso- 
nance are introduced through the following Lagrangians, 

C-^NR= ^R%AZH^d''r] + h.c., (10) 
Knr = '-^R%AYMF'^' + h.c, (11) 

^'nr= ^R%AX)d,NF'^' + h.c., (12) 

where is the field associated with the Di3(1520) resonance. Two anoma- 
lous magnetic moments k]^ and nj^ are used in conjunction with the two 
different interaction Lagrangians introduced by Benmerrouche et al. |^ for 
the 7?7i? vertex. The function O^yiV) is given by 



-^-{l + AV) + V 



V = X,Y,Z (13) 

where the parameters X, Y and Z, often referred to as off shell parameters, 
are obtained by fitting the experimental data for the elementary reaction. 

The authors of reference |Q, through fits to existing data, provided sev- 
eral sets of values for the parameters and observables related to the coupling 
constants of the effective interaction Lagrangians of Eqs. (^p!2D. The param- 
eters we use are those of the second column of table V of the above reference. 
These parameters are associated with the PS coupling which is the form used 
in the present work. The off-shell parameters a, (3 and 5 given in that table 
are related to the X, Y and Z introduced in Eqs. ( p!0| - |T^ ), as 

a = 1 + 4Z, /? = 1 + 5=l + 2X (14) 

A total of five nuclear resonances are included in the present calculations. 
Table ^ gives some of the properties of these resonances. 
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2.2 The Reaction Amplitude 

At the tree-level the S-matrix for the A( 7, rjp )A-1 reaction is |T5[ 

where \i > and |/ > are the initial and final states of the system, respectively. 
The T in front of the square bracket denotes the time ordered product of the 
operators within the bracket. We shall illustrate, below, the derivation of 
the S-matrix for one of the contributing diagrams of Fig. 2, namely the Sn 
diagram. The same procedure can be used to derive the S-matrix for the 
other diagrams. 

After introducing C^/nr and CnNB. ( from Eqs. (H), (0) and (H) ) and 
using Wick's theorem fTH], the S-matrix for the Sn resonance diagrams can 
be written as 

+V' (x) l.a^.F^" (x) R{x)R{y) ij {y) r] (y)] \i > dx^dy\ (16) 

The initial state \i > contains a photon and a target nucleus with A nucleons. 
The latter is regarded as made up of a core and a valance nucleon. We write 

|z>= al^(k^) E UJB;M,MB\Ji,Mi) 

JbMbJ 

Sj^jAJB)\"'h\„MM (17) 



where (k^) is the creation operator for a photon with momentum k^ and 

polarization s^. b^gMs creates a nucleon with angular momentum quantum 
numbers Jb and Mb bound to the core \(f)f > and S is the corresponding 
spectroscopic factor. 

The final state |/ >, composed of an rj meson, a nucleon and the residual 
nucleus, is written in a similar fashion using creation operators for the t] and 
the nucleon. To reduce the complexity of the numerical calculations we use 
the plane wave approximation for the propagators [T^. In this approximation 
the propagator for a spin-| resonance assumes the simple form 



R{x)R{y) = ^S^{x-y) = ^'§^^S\x-y), (18) 
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where S^{x — y) is the Feynman propagator of a nucleon resonance with mass 
Mn and four momentum kn. Carrying out the Fock space calculations and 
using the radiation gauge, the S-matrix of Eq. (|IBp reduces to the following 
form 



where F^^^ is a 4x4 matrix operator and contains combinations of the Dirac 
7 matrices, polarization of the interacting photon, mass and four momentum 
of the propagating resonance. The explicit form of this operator for each 
of the contributing diagrams are given below. In the above equation, E^, 
Ep and E^ are energies of the photon, proton and rj meson respectively and 
M is proton mass. \Efj^^A/g(x) describes a bound nucleon and is a solution 
of Dirac equation with the appropriate strong scalar and vector potentials. 
"^sfix) is the wavefunction of the outgoing proton which in the plane wave 
approximation is just a plane wave Dirac spinor. In the distorted wave cal- 
culations \I's^. (x) is a solution of the Dirac equation with the strong vector 
and scalar optical potentials. The rj meson wavefunction ip^ is a solution 
of the Klein-Gordon equation with an appropriate optical potential for the 
?7-nucleus system. 

The structure of the total S-matrix is the same as Sf^ of Eq. (|19D except 
for F^jj which must be replaced by a sum representing contributions from 
all the diagrams included in the model ( Fig. 2 ). We obtain the following 
operators for these diagrams 




{27rf^\Ep2Er,2Ej 

X E iJf,JB;Mf,MB\Ji,M^)\Sj^j^iJB) 



1/2 



JbMb 




(19) 



F 



proton 




(20) 




(21) 
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2 ' 



(1) 



^(1) _ UnrKr ^ 



p(2) ^ fyjNR^^ 



where 



1 „ „ 1 , 2 



and is the local momentum of the bound nucleon and 
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2.2.1 Plane Wave Calculations 



To calculate the amplitude of Eq. ( [T9| ) one can use either plane waves or 
distorted waves to describe the outgoing particles. In the plane wave ap- 
proximation the final state interactions of the outgoing nucleon and rj-meson, 
with the recoil nucleus, are ignored. The respective wavefunctions have the 
form 



E + M 
2M 

if^{x) = e-*'^''-^ " (29) 



Thus in the plane wave approximation ( PWA ) the S-matrix in Eq. (p!9|) 
takes the form 

X {Jf,JB;Mf,MB\J^,Mi)[Sj^j^.{JBf^' 
JbMb 



X 



j d^xu\f{k,)e^^^-'^Vs,,^ JbmA^)^ (30) 



where fe^ec is the momentum of the recoil nucleus ( krec = — kp — k-^i ). 
Equation (|30D shows that the PWA S-matrix is made up of terms proportional 
to the Fourier transform of the components of the bound state wavefunction. 

2.2.2 Final State Interactions: Distorted Wave Calculations 

The rest of this section is devoted to a derivation of the form of the S- 
matrix in the distorted wave approximation ( DWA ). The continuum nucleon 
wavefunction is written as 



4f{x) = Arc 



e (k/)3^Li/2j(i^) 



'L1/2J\ 

LJM 



x{L,l/2;M-Sf,Sf\J,M)[fLAr), la-vgLAr)] (31) 

where f^j and gLj refer to the upper and lower radial functions. These are 
solutions of the radial Dirac equation with vector and scalar optical potentials 
obtained from comparison with proton-nucleus elastic scattering data. More 
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details of the Dirac proton wavefunctions can be found in reference 
Expanding the eta and photon wavefunctions into partial waves and using 
these in the expression of the amplitude given by Eq. (|TP|), one gets the DWA 
S-matrix as 



e ( E„ + M\ 



1/2 



^ EpErjE^ j 



6 {Eb + E^ — Ep — E^ 



J2 iJf,JB;Mj,MB\Ji,M,)[Sj^j^.{JB] 
JbMb 

L^-L-L, ^2L^ _^ ^^1/2 



1/2 



X J2 

LJ TjYj 

E rr-' (k,) [F-" (k,) 



MMr, 



Rff [KKsLrjl^] An 
+iRfg [KKsLriLj] Ai2 

-iRgf [kkbL^L^] A21 

+ Rgg [KKBLfjLy 



K, Kb 




_M Mb 


Mr, 


Kj ^ q 




_M Mb 


M, 


' k' Kb 




_M Mb 








M Mb 


Mr, 



L,--M-Sf,Sf\J,M 



(32) 



where k = {L — J) (2 J + 1) and L' = 2J — L. The radial integrals are given 
by 

RcB [KKBLr,L^] = J Clj (r) BlbMb (r) vl„ {kr,r) Jl-, (k-yr) r'^dr, (33) 

and C (r) is either an upper- or lower- component radial function for the 
ejected proton while B (r) is the counterpart for the bound state wavefunc- 
tion. The angular integrals are given as follows 



A: 



M Mb ^] = / (^-/^^)'r.-^."?i/2..C {n)dn,i34) 
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The matrix operators Tij are the elements of the 4x4 F matrix operator 
Til 

r22 



(35) 



2.3 Observables 

It is useful to define a function of the initial and final spin projections Z^^,^^ 
by rewriting the amplitude of Eq. (^) in the following form 

1/2 



s 



fi 



5 {Ep + — Eb — E^) 



e Ep + M ' 

TT \ EpErjEy J 
JbMb 



X 



1/2 



7*/ 



(36) 



The triple differential cross section is related to -^^{/^ by 



dVl^dVlpdEr^ 



5«7r 
he 

1 



X- 



Ep + Mc^ 
Ej 



(37) 



where a is the fine structure constant and the recoil factor R is defined in a 



manner similar to the case of ( e, e'p ) reactions [|T7 



and Er is the total energy of the residual nucleus. 

The photon asymmetry for linearly polarized incident photons is 

d(7\\ — da I 

A 



da\\ + da_i 



(38) 



(39) 



where da\\ and dcr_L are the cross section for specified polarizations of the 
incident photon, namely parallel and perpendicular to the plane of the reac- 
tion. 
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3 Results and Discussion 



This paper is concerned with the study of the exclusive quasifree { J, rjp ) 
reaction on complex nuclei. To date there are no available experimental 
data for this particular reaction ( data exist only for inclusive reactions 
). Since the lack of attempts to measure the exclusive observables is mainly 
due to the expectation of low cross sections, it is of practical importance 
to identify the kinematical regions in which the cross sections are relatively 
large. Figure 3 shows the calculated cross section for this reaction on ^^C 
leading to the ground state of ^^B ( assumed to be a pure P3/2 hole state ). 
The incident photon energy is 750 MeV. The calculations are carried out in 
the plane wave approximation for the outgoing proton and 1] meson. The 
bound nucleon is described by a solution of the Dirac equation involving the 



relativistic scalar and vector Hartree potentials of reference [0. We used 



the maximum possible value for the spectroscopic factor ( i.e 2Jb + 1; Js is 
the total angular momentum of the bound proton ). The azimuthal angles 
of the outgoing proton and t] meson are fixed at 180° and 0°, respectively. 
Calculations are performed at four different polar angles for the t] meson and 
cover the outgoing proton polar angular range 10°-50°. Cross sections are 
plotted as a function of polar angle of the outgoing proton and the kinetic 
energy of the 7] meson. These results indicate that one can observe relatively 
large cross sections for 6n in the range 10°-30° and ^p=15°-30°. The pair of 
angles 9n = 9p = 30° may be a preferable setting experimentally and hence 
will be used for most of the analysis done here. Similar results were found for 
^^O and ^°Ca nuclei. These results are expected to change somewhat with 
respect to variation in the incident photon energy. The larger cross sections 
will shift towards more forward angles as the photon energy increases. 

It is well known that the eta photoproduction reaction on the proton near 
threshold is dominated by the Sii(1535) resonance formation. We have done 
calculations aimed at assessing the extent to which this dominance prevails 
when the reaction takes place within the nuclear medium. The calculations 
are done for the reaction discussed above, again in the plane wave limit, where 
the angles of the outgoing particles are fixed at On = Bp = 30°. We varied 
the energy of the incident photon from near threshold to the energy where 
the contributions from the Sii(1535) diagrams are comparable to contribu- 
tions from all other diagrams. Results of these calculations are illustrated in 
Fig. 4: the dotted curves give the cross sections due to the formation of the 



13 



Sii(1535); contributions due to the rest of the diagrams are represented by 
sohd curves. Note the varying scale for the cross section axes. Figure 4 shows 
that the Sii(1535) resonance dominates the cross section from threshold up 
to a photon energy of about 1.1 GeV. The cross section is small at ener- 
gies near threshold and increases with increasing photon energy, reaching its 
largest values in the photon energy range 750-950 MeV. At energies higher 
than 950 MeV, the cross section due to the Sii(1535) resonance decreases as 
the photon energy increases. 

In Fig. 5 we show the relative contributions of the individual diagrams. 
The cross section arising from the Sii(1535) pole diagrams, shown by the 
long-dashed curve, is clearly the dominant contribution. The second largest 
contribution comes form the proton poles( double dash-double dotted curve 
). The results of the calculations for the Di3(1520) resonance and the vector 
mesons poles are shown by dot-dot-dashed and dot- dash- dashed curves, re- 
spectively. The solid curve shows the total cross section ( i.e contributions 
from all the diagrams ). Besides confirming clearly that the Sii(1535) res- 
onance is the main contributor to the reaction in this energy region, these 
results also give us an indication of the relative importance of other diagrams. 
The next largest contribution comes from the proton pole followed by the D13 
resonance, then the vector meson diagrams. The Pii(1440), Sii(1650) and 
Pii(1710) resonances do not make significant contributions to the reaction 
at this energy. 

In order to establish that the above conclusions are not effected by the 
inclusion of the final state interactions of the outgoing particles with the 
residual nucleus, we carried out distorted wave calculations in which these 
final state interactions are taken into account. Appropriate optical potentials 
are used in the calculation of the wavefunctions of the outgoing proton and 
rj meson. In all the distorted wave calculations presented here the global 
optical potentials of Cooper et.al [|l^ have been used to calculate the Dirac 
distorted wavefunctions of the outgoing proton. For the outgoing rj meson 
we use the Klein-Gordon equation together with the optical potentials of 
references The results are shown in Fig.6; the sohd curve shows 

the cross section due to the Sii(1535) resonance, whereas the dashed curve 
includes all contributions other than those of the Sii(1535). These results 
indicate that, as in the plane wave calculations discussed above, the Sii(1535) 
resonance is the dominant contributor to the reaction in the distorted wave 
limit. Note the general suppression of the cross section as a result of the 
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inclusion of the final state interactions. 

Figure 7 shows the calculated observables for the reaction where we inves- 
tigate separately the effects of the proton and rj meson final state interactions. 



The energy dependent global optical potential of Cooper et al. |jT9[, and the 
optical potential DWl of reference |^ are used for outgoing proton and r] 
meson, respectively. Plane wave and various distorted wave calculations are 
shown. The comparisons of the cross sections are shown in Fig. 7- (a) and 
those for the photon asymmetry are shown in Fig. 7-(b). The short dashed 
curve shows the suppression due to the final state interaction of the outgoing 
proton with the residual nucleus. This suppression is more evident at higher 
7] energies. These correspond to smaller energies of the outgoing protons and 
indicate that these protons are more affected by the distortion. This same 
qualitative behavior holds for the f] meson distortion ( dashed curve ); the 
suppression is skewed towards smaller r] energies. The dotted curve is the 
calculated cross section when distortion effects are included for both of the 
outgoing particles. The suppression of the cross section in this case is much 
more pronounced. 

Figure 7-(b) shows the photon asymmetry for the same reaction ( see 
Eq. (|3^ )). Plane wave calculations (solid curve) produce a nearly fiat curve 
for the photon asymmetry. Distortion of the outgoing rj meson affects the 
shape and magnitude of the photon asymmetry only slightly. Inclusion of the 
final state interactions of the outgoing proton results in larger changes in the 
asymmetry. When both final state interactions are included, the resulting 
curve ( dotted curve ) has characteristics close to the case when only the 
proton is distorted. 

Figure 8 shows the effects of different rj optical potentials on the calculated 
observables. The curves of Fig. 8 are obtained with the same bound and 
continuum proton potentials while different optical potentials are used for 
the outgoing 77 meson. Five different potentials are used. The first two are 
those (denoted by DWl and DW2) used by Lee et al. [|]. These are obtained 
in a tp approximation using rjN t-matrices obtained in a coupled channel 
approach [7|, ^ ( a recent work by Batinic et al. |2T[] points out to a sign 
error in their analysis in which may mean that the potential DW2 is 
not reliable ). The other sets of potentials are those of Chiang et al. [p^] . 
These are obtained by calculating the 77 self energy in nuclear matter taking 
into account only contributions from the Sn resonance. The three different 
potentials used correspond to three different choices of the real part of the 
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N* self energy ( namely ReSAr* = {p/po)Vn* with VAr.=50, 0, -50 MeV; p and 
Po are the nuclear density for finite and infinite nuclear matter, respectively. 
) The uncertainty about the parameter V^* is a drawback for the potentials 
of Chiang et ai. 

In general the calculated cross sections show some dependence on the 
type of distorting potential used; there are some variations in shape and 
magnitude. The magnitude of the cross section obtained with the DWl 
potential is somewhat larger than those predicted by the other potentials 
near the mid-range of i] energies. 

The corresponding calculations for photon asymmetry with different opti- 
cal potentials for the outgoing t] meson are shown in Fig.8-(b). All potentials 
produce approximately the same shape and magnitude for the photon asym- 
metry. The asymmetry therefore, at least at this energy, is not particularly 
sensitive to the different choices of the final state interactions of 77 mesons. 

We have also calculated the observables of the reaction using several dif- 
ferent proton optical potentials. We found that the calculated observables 
are rather insensitive to the variations among these optical potentials. 

The sensitivity of the results to the bound state potentials was also as- 
sessed by performing calculations in which different binding potentials are 
used. In this case relativistic scalar and vector potentials of Woods Saxon 
shape 1^ were used to generate the Dirac bound state wavefunction. The 
results are shown by the dashed curve in Fig.9-(a) along with the results 
obtained using the Hartree potentials adopted throughout the present study 
( solid curve ). Comparison of these curves indicates the level of sensitivity 
in our model calculations to the bound state potentials. This comparison 
shows that the Hartree bound state potentials lead to somewhat larger cross 
sections at this photon energy. On the other hand the photon asymmetry 
calculations appear to be insensitive to different choices of the bound state 
potentials (Fig.9-(b)). 

In addition to the kinematics used above, we also calculated the observ- 
ables of the reaction for another set of angles for the outgoing particles. The 
motivation for choosing this new set is to compare the results of our model 
with those of the nonrelativistic model of Lee et al. [§. Figure 10 shows this 
comparison. The calculations are performed at the following angular settings 
for the outgoing particles: = 20°, 0^ = 0° and 6^ = 15°, 0p = 180°. We 
have plotted the results of the non-relativistic plane wave calculations of Lee 
et al. M together with our results for both plane wave and distorted wave 
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calculations. Figure 10- (a) shows the cross section and 10- (b) the result- 
ing photon asymmetry. Results of the relativistic plane wave and distorted 
wave calculations are shown by solid and short dashed curves, respectively. 
The nonrelativistic plane wave calculations are depicted by the long dashed 
curves. 

Comparison of the solid and long dashed curves indicates that both rel- 
ativistic and nonrelativistic calculations produce curves with similar shape, 
but the nonrelativistic calculations predict generally larger cross sections in 
the f] energy range up to ~120 MeV, and fall below at higher rj energies. 
There is also a slight shift in the relative positions of the peaks in the two 
calculations. As noted earlier, the inclusion of final state interactions leads 
to a suppression of the calculated cross sections. 

The plane wave relativistic and nonrelativistic predictions for the photon 
asymmetry differ significantly ( Note that Lee et al. use a definition of 
asymmetry which is opposite in sign to the one used here; for the purpose 
of the present comparison the results from reference have been multiplied 
by -1 ). The relativistic calculations produce a flat curve with small positive 
values, whereas the nonrelativistic calculations result in a curve with large 
negative values, with strong dependence on the rj energy. The relativistic 
calculations show some sensitivity to final state interactions whereas Lee et 
al. report such effects are negligible in the nonrelativistic calculations. 

We carried out similar calculations for the photoproduction reaction on a 
^^O target leading to the ground state in ^^N. Figure 11 shows the results of 
these calculations. The residual nuclear state is assumed to be a pure lpi/2 
hole state. The Hartree potentials of reference are used to calculate 



the bound state wavefunction. The final state interactions of the outgoing 
proton and rj meson are calculated using the energy dependent global optical 



potential of reference [^] and the DWl optical potential of reference [g], 
respectively. 

The results for the cross section are very much similar to the ^^C case 
( see Fig. 7 ). In the case of ^^O we note that the distorted wave cross section 
( dotted curve ) is slightly smaller in comparison with the corresponding case 
in ^^C ( dotted curve in Fig. 7 ). On the other hand the resulting photon 
asymmetries for this target, behave differently form those of the ^^C target. 
The distorted wave calculations here peak at lower rj energies. Moreover 
asymmetries for this target have larger magnitudes but this occurs in a region 
of very small cross sections. These differences with the ^^C case indicate 
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that the asymmetry is dependent on the single particle state from which the 
proton is ejected. 

We have also used our relativistic model for calculations of the reaction 
on ^''Ca. Figure 12 shows the results of our calculations for the case in which 
the residual nucleus ( ^^K ) is left in its ground state ( a ld3/2 hole state ) 
and the incident photon energy is 750 MeV. The calculated cross sections of 
Fig. 12- (a) show curves with different shapes from those of the two targets 
studied previously. The final state interactions of the outgoing particles lead 
to a suppression of the cross section for this target in a manner similar to 
what has been noted earlier for ^^C and ^^O. We find that the earlier obser- 
vation on photon asymmetries is confirmed in the present case; the photon 
asymmetry calculations are sensitive to the final state interactions of the out- 
going particles. This is in contrast to the predictions of the nonrelativistic 
model for the reaction |Q . The shape and magnitude of the asymmetries are 
different from those found for the ^^C and ^^O targets. 

We have also carried out calculations for the observables of the reaction 
on a ^°Ca target leading to an odd parity lp3/2 hole state in ^^K. The results 
are shown in Fig. 13. The curves in this figure have the same features as in 
Fig. 10, except that the differences in the magnitude of the relativistic and 
nonrelativistic cross sections are larger. Also the effects of the final state 
interactions are more prominent for the present case. Note that the photon 
asymmetries given in Figs. lO-(b) and 13-(b) have approximately the same 
shapes. These similarities and the fact that reaction is taking place on a 
proton in a lp3/2 state, support our earlier statement on the state dependency 
of the asymmetry. 

The calculations reported in this section have all used a specific set of 
coupling constants; one of several that are given in reference Q . The question 
then arises concerning the sensitivity of the calculated results to the possible 
differences among these sets which provide nearly equivalent fits to the data 
for the elementary reactions. We have carried out calculations using the set 
given in the first column of table V of reference The results show that 
there are some changes in the magnitudes of the calculated cross sections 
and photon asymmetries. These changes are at the level of 15% for the cross 
sections but are somewhat larger for the asymmetries. 
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4 Conclusions 



In this paper we have developed a relativistic model for the quasifree photo- 
production of 7] mesons on complex nuclei. The ingredients of the model are: 
i) the nucleon wave functions are solutions of the Dirac equation with appro- 
priate scalar and vector potentials consistent with the relativistic mean field 
approach, ii) the rj meson is described by solutions of the Klein-Gordon equa- 
tion with appropriate optical potentials, and iii) the interactions between the 
fields are introduced through a covariant effective Lagrangian. Contributions 
from nucleon resonances, t-channel vector meson as well as the nucleon Born 
diagrams are included. The model is used to calculate the cross section and 
photon asymmetry, for different target nuclei. Both plane wave and distorted 
wave calculations are presented. 

We wish to comment on the question of the gauge invariance of the reac- 
tion amplitude calculated in the above manner. The effective Lagrangian we 
start with is gauge invariant. As we look at the nuclear amplitude itself, as 
given by Eq. ([T9|) , we find that its invariance is determined by the behavior 
of the r operators given in Eqs.(^D[]2^). All these operators are gauge invari- 
ant, except for the electric parts of Eq. (]20|) . These latter terms are invariant 
only in the limit where the particles involved in the elementary interaction 
are considered to be on their mass shell. 

The energy region where the Sii(1535) resonance dominates the reaction 
is identified. As expected this region is close to threshold and covers a range 
of photon energies from 750 MeV to 950 MeV. None of the other diagrams 
makes contributions at the same level as the Sn. We find that the next 
leading diagrams are proton, D13 resonance and vector meson poles. However 
we also find that there is some cancellation among these contributions. In our 
view the energy region identified above is perhaps the best suited for using 
the rj photoproduction reactions to study the properties of the Sn resonance 
in the nuclear medium. 

When the final state interactions of the outgoing particles are included, 
we find that the resulting cross sections are strongly suppressed. The extent 
of this suppression depends somewhat on the type of optical potential used 
to represent the final state interaction. We find some variations among the 
available sets of potentials particularly those that describe the r] meson final 
state interactions. Several r] optical potentials were tested and lead to dif- 
ferent results for the cross section as well as the photon asymmetry . The 
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uncertainty in rj optical potentials makes it difficult to make final predictions 
for the cross sections. More theoretical work on the final state interactions 
of the r] is required. 

Calculations were carried out using the same kinematics as those in the 
recent nonrelativistic calculations of Lee et. al |@] . Results of our calculations 
for the cross section have shapes close to those of these nonrelativistic calcu- 
lations but the magnitudes of our results are somewhat smaller. In contrast, 
the photon asymmetries predicted by the two models differ significantly. The 
nonrelativistic model gives large asymmetries in the plane wave limit, which 
are found to be insensitive to the final state interactions of the outgoing parti- 
cles. Our predictions yield much smaller asymmetries and these are strongly 
affected by the final state interactions. Lee et al. have also found that a 
small percentage change in the masses of the Sn and Di^ have large effects 
on the cross sections and the asymmetries, respectively. Our present calcu- 
lations show similar effects. These are indications that eta photoproduction 
reactions can be used as a probe of medium modifications of the properties 
of these resonances. These of course are expected to go beyond the simple 
mass shifts; for example there may in fact be a radial density dependence 
of the masses of the propagators. Studies of these effects with the present 
relativistic model are planned. 

The above statements indicate that the present reaction is rich in its 
physics content. If enough data become available, one will be able to clarify 
the role played by the intermediate nucleon states ( the Sn and the D13 ) 
as well as the vector mesons while embedded in the nuclear environment. 
The question of the optical potential for the 77-nucleus system could also 
be clarified. The data could also be used to assess the need for a strictly 
relativistic treatment. It is clear that there is a multiplicity of effects at 
play and hence it is imperative that measurements at several energies in 
the region discussed above and for a variety of target nuclei be carried out. 
Measurements of the cross sections and photon asymmetries in exclusive 
production, as well as coherent production, will be very helpful in shedding 
more light on the above issues. In a subsequent paper the present model will 
be extended to the treatment of inclusive reactions for which some data have 
become available recently [Q. 
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Table Caption 



TABLE 1. Summery of the properties of the baryon resonances included 
in present model. J'^ is the spin parity, L is the relative angular momentum 
related to ttN scattering state. The subscripts I and J for the L values ( third 
column ) represent the isospin and total angular momentum assigned to the 
resonance, respectively. F is the total width of the resonance. 



Figure Captions 

FIG. 1. Feynman diagram for an A{ 7, rjp )A — 1 reaction. The incident 
photon with momentum k is absorbed by a proton embedded in nucleus A, 
resulting in ejection of a proton and an rj meson. 

FIG. 2. The contributing Feynman diagrams to A{ 7, rip)A — l reaction, a ) 
and b) are the s- and u-channel Born diagram respectively, c) the t-channel 
vector meson diagram, d) and e) are the s- and u-channel nucleon resonance 
poles. 

FIG. 3. The cross section for the ^^C( 7, r/p Y^Bg ^^ reaction for the regions 
of the phase space where the reaction has significant yield. Calculations are 
performed in the plane wave approximation. 

FIG. 4. The plane wave cross sections for the ^^C( 7, //p Y^Bg^ reaction for 
photon energies from near threshold up to 1.2 GeV. 

FIG. 5. Contributions of different diagrams to the cross section of the 
^^C( 7,77P Y'^Bgs reaction at photon energies of 750 MeV. 

FIG. 6. Differential cross section for the same reaction as Fig. 5, calculated 
in the distorted wave approximation. Solid curve - DW calculations using 
only Sii(1535) diagrams ( curve labelled Sn ). Long dashed curve - DW 
calculations using all the diagrams but those of Sii(1535) resonance ( curve 
labelled Rest ). 

FIG. 7. Differential cross section (a), and photon asymmetry (b) for the reac- 
tion ^^C(7, rjpY^Bg s. at = 750 MeV. The Hartree potentials of reference 



ISl are used in calculation of the bound state wavefunction. The final state 
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energy dependent global optical potentials are taken from reference ||19||. The 
7] optical potential is the DWl potential of Lee et al. [§]. Solid curve - plane 
wave calculations. Long dashed curve - calculations include only final state 
interactions of the rj meson ( r] Distorted ) . Short dashed curve - calculations 
include only final state interactions of outgoing proton ( Proton Distorted ). 
Dotted curve - both rj and proton waves are distorted ( Full DW ). 

FIG. 8. Same as Fig. 7 but for different choices of rj optical potentials. Solid 
curve - DW calculations using DWl of reference . Long dashed curve - DW 
calculations using DW2 of reference [Q. DW calculations using the optical 
potential of reference with the real part of the Sn self energy set to 50 
MeV (short dashed curve), Mev ( dotted curve) and -50 MeV ( dash-dotted 
curve). 

FIG. 9. Differential cross section (a), and photon asymmetry (b) for the 
same reaction as in Fig. 7. The final state optical potentials for the rj meson 
and outgoing proton are the same as those of Fig. 7. Curves are labelled for 
different choices of proton bound potentials: Solid curve - DW calculations 
using Hartree potential of reference [|18| and dashed curve - DW calculations 



using Woods-Saxon potential of reference . 

FIG. 10. Differential cross section (a), and photon asymmetry (b) for the 
same reaction as Fig. 7. Curves are labelled as: Solid curve - relativistic PW 
calculations. Long dashed curve - nonrelativistic PW calculations of Lee et 
al. Short dashed curve - relativistic DW calculations. 

FIG. 11. Differential cross section (a), and photon asymmetry (b) for the 
reaction ^^0(7, f]pY^Ng s_ at E.y = 750 MeV. The Hartree potential of refer- 
ence HTBH is used in calculation of the bound state wavefunction. The final 



state energy dependent global optical potentials are taken from reference |T9 
The r] optical potential is the DWl potential of Lee et al. Solid curve - 
plane wave calculations. Long dashed curve - calculations include only final 
state interactions of the rj meson with nuclei ( rj Distorted ). Short dashed 
curve - calculations include only final state interactions of outgoing proton 
( Proton Distorted ). Dotted curve - both r] and proton waves are distorted 
( DW ). 

FIG. 12. Differential cross section (a), and photon asymmetry (b) for the 
reaction '^^Ca{'j, TjpY^Kg s. at = 750 MeV. The Hartree potential of refer- 
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ence [|T8l is used in the calculation of the bound state wavefunction. The final 



state energy dependent global optical potentials are taken from reference ||T9 
The rj optical potential is the DWl potential of Lee et al. 0. Solid curve - 
plane wave calculations. Long dashed curve - calculations include only final 
state interactions of the rj meson with nuclei ( r] Distorted ). Short dashed 
curve - calculations include only final state interactions of outgoing proton 
( Proton Distorted ). Dotted curve - both rj and proton waves are distorted 
( DW ). 

FIG. 13. Differential cross section (a), and photon asymmetry (b) for the 
same reaction as Fig. 12. The reaction is assumed to take place in lp3/2 
proton. The potentials are the same as those of Fig. 12. Curves are labelled 
as in Fig. 10. 
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Resonance 


J'^ 


L2/2J 


Mass(MeV) 


r(MeV) 


N*(1440) 


1/2+ 


Pii 


1440 


350 


N*(1520) 


3/2" 


Di3 


1520 


120 


N*(1535) 


1/2" 


Sii 


1535 


150 


N*(1650) 


1/2- 


Sii 


1650 


150 


N*(1710) 


1/2+ 


Pii 


1710 


100 



Table 1: 
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